-------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\JOP IRT Brazil Res
> ults.log
  log type:  text
 opened on:   6 Mar 2018, 12:11:18

. use "`stick'\Data\Brazil\BEPS 2014\BEPS 2014 Merged Data v3.dta", clear

. *use "C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\Data\Brazil\BEPS 2014\BEPS 2014 Merged Data v3.dta", clear
. 
. xtdes

idtelefone:  1, 2, ..., 3120                                 n =       3120
    wave:  1, 2, ..., 7                                      T =          7
           Delta(wave) = 1 unit
           Span(wave)  = 7 periods
           (idtelefone*wave uniquely identifies each observation)

Distribution of T_i:   min      5%     25%       50%       75%     95%     max
                         1       1       1         2         4       5       5

     Freq.  Percent    Cum. |  Pattern
 ---------------------------+---------
     1355     43.43   43.43 |  1......
      149      4.78   48.21 |  1.1.111
      137      4.39   52.60 |  11..111
      131      4.20   56.79 |  1.11.11
      120      3.85   60.64 |  11.1.11
       92      2.95   63.59 |  11.....
       87      2.79   66.38 |  1....1.
       80      2.56   68.94 |  1..1.11
       74      2.37   71.31 |  1...111
      895     28.69  100.00 | (other patterns)
 ---------------------------+---------
     3120    100.00         |  XXXXXXX

. 
. fillin idtelefone wave

. sort idtelefone wave

. 
. recode _fillin (1=0) (0=1), gen(respond)
(21840 differences between _fillin and respond)

. 
. 
. *I. PARTY ID MEASURES
.         *1: Working with the "sympathy" variable that was asked of everyone in every wave
. recode vb10 (.a=.) (.b=.) (.c=.)
(vb10: 84 changes made)

. recode vb11 (.a=.) (.b=.) (.c=.)
(vb11: 5849 changes made)

. replace vb11=. if vb10==.
(1 real change made, 1 to missing)

. 
. gen symp=vb11
(19,885 missing values generated)

. replace symp=0 if vb10==2
(5,725 real changes made)

. replace symp=. if vb10==.
(0 real changes made)

. 
. xi, noomit: ci proportions i.symp if wave==1

                                                         -- Binomial Exact --
    Variable |        Obs  Proportion    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
    _Isymp_0 |      3,042    .7945431    .0073255        .7797408    .8087766
 _Isymp_1501 |      3,042    .1180145    .0058495        .1067567    .1300125
 _Isymp_1502 |      3,042     .026956    .0029364         .021495    .0333502
 _Isymp_1503 |      3,042    .0259698    .0028836        .0206132    .0322621
 _Isymp_1504 |      3,042    .0075608    .0015706        .0047988    .0113235
-------------+---------------------------------------------------------------
 _Isymp_1505 |      3,042    .0016437    .0007345        .0005339    .0038315
 _Isymp_1506 |      3,042    .0016437    .0007345        .0005339    .0038315
 _Isymp_1507 |      3,042    .0039448    .0011365        .0020399    .0068806
 _Isymp_1508 |      3,042    .0062459    .0014284        .0037645    .0097366
 _Isymp_1509 |      3,042     .003616    .0010883        .0018065    .0064609
-------------+---------------------------------------------------------------
 _Isymp_1510 |      3,042    .0098619    .0017916        .0066634    .0140488

. 
. recode vb10 (2=1) (1=0) (.=.), gen(symp_indep)
(7721 differences between vb10 and symp_indep)

. 
. recode vb11 (1501=1) (else=0), gen(symp_pt)
(21840 differences between vb11 and symp_pt)

. recode vb11 (1502=1) (else=0), gen(symp_pmdb)
(21840 differences between vb11 and symp_pmdb)

. recode vb11 (1503=1) (else=0), gen(symp_psdb)
(21840 differences between vb11 and symp_psdb)

. recode vb11 (1504=1) (else=0), gen(symp_psb)
(21840 differences between vb11 and symp_psb)

. recode vb11 (1505=1) (else=0), gen(symp_dem)
(21840 differences between vb11 and symp_dem)

. recode vb11 (1506=1) (else=0), gen(symp_ptb)
(21840 differences between vb11 and symp_ptb)

. recode vb11 (1507=1) (else=0), gen(symp_psol)
(21840 differences between vb11 and symp_psol)

. recode vb11 (1508=1) (else=0), gen(symp_pv)
(21840 differences between vb11 and symp_pv)

. recode vb11 (1509=1) (else=0), gen(symp_pdt)
(21840 differences between vb11 and symp_pdt)

. 
. replace symp_indep=. if vb10==. 
(0 real changes made)

. replace symp_pt=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_pmdb=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_psdb=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_psb=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_dem=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_ptb=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_psol=. if vb10==.
(14,119 real changes made, 14,119 to missing)

. replace symp_pv=. if vb10==. 
(14,119 real changes made, 14,119 to missing)

. replace symp_pdt=. if vb10==.
(14,119 real changes made, 14,119 to missing)

. 
. replace symp_indep=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_pt=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_pmdb=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_psdb=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_psb=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_dem=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_ptb=. if vb11==. & vb10==1 
(41 real changes made, 41 to missing)

. replace symp_psol=. if vb11==. & vb10==1
(41 real changes made, 41 to missing)

. replace symp_pv=. if vb11==.  & vb10==1
(41 real changes made, 41 to missing)

. replace symp_pdt=. if vb11==. & vb10==1
(41 real changes made, 41 to missing)

. 
. label define pid 0 "Independent" 1501 "PT" 1502 "PMDB" 1503 "PSDB" 1504 "PSB" 1505 "DEM" 1506 "PTB" 1507 "PSOL" 1508 "PV" 1509 "PDT" 1510 "Other"

. label values symp pid

. 
. tab1 symp* if wave==1

-> tabulation of symp if wave==1 

       symp |      Freq.     Percent        Cum.
------------+-----------------------------------
Independent |      2,417       79.45       79.45
         PT |        359       11.80       91.26
       PMDB |         82        2.70       93.95
       PSDB |         79        2.60       96.55
        PSB |         23        0.76       97.30
        DEM |          5        0.16       97.47
        PTB |          5        0.16       97.63
       PSOL |         12        0.39       98.03
         PV |         19        0.62       98.65
        PDT |         11        0.36       99.01
      Other |         30        0.99      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_indep if wave==1 

  RECODE of |
vb10 (vb10. |
 atualmente |
o(a) sr.(a) |
  simpatiza |
  com algum |
    partido |
    politic |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        625       20.55       20.55
          1 |      2,417       79.45      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_pt if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      2,683       88.20       88.20
          1 |        359       11.80      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_pmdb if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      2,960       97.30       97.30
          1 |         82        2.70      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_psdb if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      2,963       97.40       97.40
          1 |         79        2.60      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_psb if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,019       99.24       99.24
          1 |         23        0.76      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_dem if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,037       99.84       99.84
          1 |          5        0.16      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_ptb if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,037       99.84       99.84
          1 |          5        0.16      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_psol if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,030       99.61       99.61
          1 |         12        0.39      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_pv if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,023       99.38       99.38
          1 |         19        0.62      100.00
------------+-----------------------------------
      Total |      3,042      100.00

-> tabulation of symp_pdt if wave==1 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      3,031       99.64       99.64
          1 |         11        0.36      100.00
------------+-----------------------------------
      Total |      3,042      100.00

. tab1 symp* if wave==6

-> tabulation of symp if wave==6 

       symp |      Freq.     Percent        Cum.
------------+-----------------------------------
Independent |        884       73.48       73.48
         PT |        200       16.63       90.11
       PMDB |         29        2.41       92.52
       PSDB |         42        3.49       96.01
        PSB |         13        1.08       97.09
        DEM |          1        0.08       97.17
        PTB |          3        0.25       97.42
       PSOL |          5        0.42       97.84
         PV |         11        0.91       98.75
        PDT |          5        0.42       99.17
      Other |         10        0.83      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_indep if wave==6 

  RECODE of |
vb10 (vb10. |
 atualmente |
o(a) sr.(a) |
  simpatiza |
  com algum |
    partido |
    politic |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        319       26.52       26.52
          1 |        884       73.48      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_pt if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,003       83.37       83.37
          1 |        200       16.63      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_pmdb if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,174       97.59       97.59
          1 |         29        2.41      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_psdb if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,161       96.51       96.51
          1 |         42        3.49      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_psb if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,190       98.92       98.92
          1 |         13        1.08      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_dem if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,202       99.92       99.92
          1 |          1        0.08      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_ptb if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,200       99.75       99.75
          1 |          3        0.25      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_psol if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,198       99.58       99.58
          1 |          5        0.42      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_pv if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,192       99.09       99.09
          1 |         11        0.91      100.00
------------+-----------------------------------
      Total |      1,203      100.00

-> tabulation of symp_pdt if wave==6 

  RECODE of |
vb11 (vb11. |
   com qual |
    partido |
       o(a) |
    sr.(a). |
 simpatiza? |
          ) |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,198       99.58       99.58
          1 |          5        0.42      100.00
------------+-----------------------------------
      Total |      1,203      100.00

. 
. recode symp_indep (1=0) (0=1), gen(symp_partisan)
(7680 differences between symp_indep and symp_partisan)

. 
.         *2: Working with the "preference" variable, treatment A
. rename vb10expa exp_1

. recode exp_1 (.a=.) (.b=.) (.c=.) (77=0), gen(preference)
(3759 differences between exp_1 and preference)

. tab exp_1 if wave==1

       exp 1. qual e o seu |
       partido politico de |
              preferencia? |      Freq.     Percent        Cum.
---------------------------+-----------------------------------
            nenhum/nao tem |        562       55.70       55.70
                pt/petista |        278       27.55       83.25
          pmdb/pemedebista |         51        5.05       88.31
  psdb/pessedebista/tucano |         59        5.85       94.15
            psb/pessebista |         12        1.19       95.34
               dem/demista |          5        0.50       95.84
             ptb/petebista |          5        0.50       96.33
            psol/pesolista |          7        0.69       97.03
                pv/pevista |         12        1.19       98.22
             pdt/pedetista |          7        0.69       98.91
outro, nomes e referencias |         11        1.09      100.00
---------------------------+-----------------------------------
                     Total |      1,009      100.00

. tab symp preference

            |                          RECODE of exp_1 (exp 1. qual e o seu partido politico de preferencia?)
       symp |         0       1501       1502       1503       1504       1505       1506       1507       1508       1509       1510 |     Total
------------+-------------------------------------------------------------------------------------------------------------------------+----------
Independent |       819        150         23         29          8          1          4          5          7          5          7 |     1,058 
         PT |         8        190          2          1          1          2          0          0          1          0          0 |       205 
       PMDB |         1          1         36          1          0          0          0          0          0          0          0 |        39 
       PSDB |         1          1          0         41          1          0          1          0          0          0          0 |        45 
        PSB |         0          1          0          0          5          0          0          0          0          0          0 |         6 
        DEM |         0          0          0          0          0          2          0          0          0          0          1 |         3 
        PTB |         0          0          1          0          0          0          3          0          0          0          0 |         4 
       PSOL |         1          0          0          0          0          0          0          5          0          0          0 |         6 
         PV |         1          0          0          0          0          0          0          1         10          0          0 |        12 
        PDT |         0          1          1          0          0          0          0          0          0          3          0 |         5 
      Other |         4          3          0          1          0          0          0          0          0          0          5 |        13 
------------+-------------------------------------------------------------------------------------------------------------------------+----------
      Total |       835        347         63         73         15          5          8         11         18          8         13 |     1,396 


. label values preference pid

. 
. xi, noomit: ci proportions i.preference if wave==1

                                                         -- Binomial Exact --
    Variable |        Obs  Proportion    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
_Iprefere~_0 |      1,009    .5569871    .0156381        .5257126    .5879282
_Iprefe~1501 |      1,009    .2755203    .0140651        .2481458    .3042112
_Iprefe~1502 |      1,009    .0505451    .0068965        .0378624     .065925
_Iprefe~1503 |      1,009    .0584737    .0073867        .0448074    .0747819
_Iprefe~1504 |      1,009     .011893    .0034127          .00616    .0206824
-------------+---------------------------------------------------------------
_Iprefe~1505 |      1,009    .0049554    .0022106        .0016109    .0115261
_Iprefe~1506 |      1,009    .0049554    .0022106        .0016109    .0115261
_Iprefe~1507 |      1,009    .0069376     .002613        .0027937    .0142415
_Iprefe~1508 |      1,009     .011893    .0034127          .00616    .0206824
_Iprefe~1509 |      1,009    .0069376     .002613        .0027937    .0142415
-------------+---------------------------------------------------------------
_Iprefe~1510 |      1,009    .0109019    .0032691        .0054544    .0194226

. 
. recode preference (0=1) (else=0), gen(pref_indep)
(21840 differences between preference and pref_indep)

. recode preference (1501=1) (else=0), gen(pref_pt)
(20990 differences between preference and pref_pt)

. recode preference (1502=1) (else=0), gen(pref_pmdb)
(20990 differences between preference and pref_pmdb)

. recode preference (1503=1) (else=0), gen(pref_psdb)
(20990 differences between preference and pref_psdb)

. recode preference (1504=1) (else=0), gen(pref_psb)
(20990 differences between preference and pref_psb)

. recode preference (1505=1) (else=0), gen(pref_dem)
(20990 differences between preference and pref_dem)

. recode preference (1506=1) (else=0), gen(pref_ptb)
(20990 differences between preference and pref_ptb)

. recode preference (1507=1) (else=0), gen(pref_psol)
(20990 differences between preference and pref_psol)

. recode preference (1508=1) (else=0), gen(pref_pv)
(20990 differences between preference and pref_pv)

. recode preference (1509=1) (else=0), gen(pref_pdt)
(20990 differences between preference and pref_pdt)

. 
. replace pref_indep=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_pt=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_pmdb=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_psdb=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_psb=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_dem=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_ptb=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_psol=. if preference==.
(20,422 real changes made, 20,422 to missing)

. replace pref_pv=. if preference==. 
(20,422 real changes made, 20,422 to missing)

. replace pref_pdt=. if preference==.
(20,422 real changes made, 20,422 to missing)

. 
. recode pref_indep (1=0) (0=1), gen(pref_partisan)
(1418 differences between pref_indep and pref_partisan)

. 
.         *3: Working with the "identity" variable, treatment B
. rename vb10expb exp_2

. recode exp_2 (.a=.) (.b=.) (.c=.) (77=0), gen(identity)
(3897 differences between exp_2 and identity)

. tab identity if wave==1

  RECODE of |
 exp_2 (exp |
 2. voce se |
  considera |
 pedetista, |
pemedebista |
          , |
pessebista, |
         pe |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        660       67.28       67.28
       1501 |        180       18.35       85.63
       1502 |         35        3.57       89.19
       1503 |         46        4.69       93.88
       1504 |         20        2.04       95.92
       1505 |          3        0.31       96.23
       1506 |          3        0.31       96.53
       1507 |          5        0.51       97.04
       1508 |          9        0.92       97.96
       1509 |          5        0.51       98.47
       1510 |         15        1.53      100.00
------------+-----------------------------------
      Total |        981      100.00

. tab symp identity

            |                     RECODE of exp_2 (exp 2. voce se considera pedetista, pemedebista, pessebista, pe
       symp |         0       1501       1502       1503       1504       1505       1506       1507       1508       1509       1510 |     Total
------------+-------------------------------------------------------------------------------------------------------------------------+----------
Independent |       900         96         20         22         11          3          2          4          4          1          4 |     1,067 
         PT |        35        143          0          2          0          0          0          0          0          1          1 |       182 
       PMDB |         6          4         29          1          0          0          0          0          0          0          0 |        40 
       PSDB |         9          1          0         26          0          0          0          1          0          0          0 |        37 
        PSB |         1          0          0          1         13          0          0          0          0          0          0 |        15 
        PTB |         0          0          0          0          0          0          1          0          0          0          0 |         1 
       PSOL |         0          0          0          1          0          0          0          5          0          0          0 |         6 
         PV |         4          0          0          0          0          0          0          0          6          0          0 |        10 
        PDT |         2          0          0          0          1          0          0          0          0          5          0 |         8 
      Other |         2          1          0          0          0          0          0          0          0          0         13 |        16 
------------+-------------------------------------------------------------------------------------------------------------------------+----------
      Total |       959        245         49         53         25          3          3         10         10          7         18 |     1,382 


. label values identity pid

. 
. xi, noomit: ci proportions i.identity if wave==1

                                                         -- Binomial Exact --
    Variable |        Obs  Proportion    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
_Iidentity_0 |        981    .6727829    .0149803        .6424267    .7020977
_Iident~1501 |        981    .1834862     .012358        .1597341    .2091547
_Iident~1502 |        981    .0356779    .0059221        .0249744    .0492721
_Iident~1503 |        981    .0468909    .0067496        .0345314    .0620533
_Iident~1504 |        981    .0203874     .004512        .0124965    .0313117
-------------+---------------------------------------------------------------
_Iident~1505 |        981    .0030581    .0017629        .0006311    .0089108
_Iident~1506 |        981    .0030581    .0017629        .0006311    .0089108
_Iident~1507 |        981    .0050968    .0022736        .0016569    .0118539
_Iident~1508 |        981    .0091743     .003044        .0042034     .017344
_Iident~1509 |        981    .0050968    .0022736        .0016569    .0118539
-------------+---------------------------------------------------------------
_Iident~1510 |        981    .0152905    .0039177        .0085826    .0250942

. 
. recode identity (0=1) (else=0), gen(iden_indep)
(21840 differences between identity and iden_indep)

. recode identity (1501=1) (else=0) , gen(iden_pt)
(20861 differences between identity and iden_pt)

. recode identity (1502=1) (else=0), gen(iden_pmdb)
(20861 differences between identity and iden_pmdb)

. recode identity (1503=1) (else=0), gen(iden_psdb)
(20861 differences between identity and iden_psdb)

. recode identity (1504=1) (else=0), gen(iden_psb)
(20861 differences between identity and iden_psb)

. recode identity (1505=1) (else=0), gen(iden_dem)
(20861 differences between identity and iden_dem)

. recode identity (1506=1) (else=0), gen(iden_ptb)
(20861 differences between identity and iden_ptb)

. recode identity (1507=1) (else=0), gen(iden_psol)
(20861 differences between identity and iden_psol)

. recode identity (1508=1) (else=0), gen(iden_pv)
(20861 differences between identity and iden_pv)

. recode identity (1509=1) (else=0), gen(iden_pdt)
(20861 differences between identity and iden_pdt)

. 
. replace iden_indep=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_pt=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_pmdb=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_psdb=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_psb=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_dem=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_ptb=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_psol=. if identity==.
(20,431 real changes made, 20,431 to missing)

. replace iden_pv=. if identity==. 
(20,431 real changes made, 20,431 to missing)

. replace iden_pdt=. if identity==.
(20,431 real changes made, 20,431 to missing)

. 
. recode iden_indep (1=0) (0=1), gen(iden_partisan)
(1409 differences between iden_indep and iden_partisan)

. 
.         *4: Working with the "way of thinking" variable, treatment C
. rename vb10expc1 exp_3a

. rename vb10expc2 exp_3b

. recode exp_3a (.a=.) (.b=.) (.c=.)
(exp_3a: 3061 changes made)

. recode exp_3b (.a=.) (.b=.) (.c=.)
(exp_3b: 4038 changes made)

. gen thinking=exp_3b
(21,551 missing values generated)

. replace thinking=0 if exp_3a==2
(973 real changes made)

. 
. label values thinking pid

. 
. xi, noomit: ci proportions i.thinking if wave==1

                                                         -- Binomial Exact --
    Variable |        Obs  Proportion    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
_Ithinking_0 |        918    .7908497    .0134232        .7630826     .816736
_Ithink~1501 |        918    .1045752    .0100997        .0855294    .1261998
_Ithink~1502 |        918    .0250545    .0051584        .0159472    .0373579
_Ithink~1503 |        918    .0217865    .0048182        .0133573    .0334477
_Ithink~1504 |        918    .0098039    .0032519        .0044925     .018529
-------------+---------------------------------------------------------------
_Ithink~1505 |        918     .003268    .0018837        .0006744    .0095204
_Ithink~1506 |        918    .0043573    .0021739        .0011885    .0111185
_Ithink~1507 |        918    .0043573    .0021739        .0011885    .0111185
_Ithink~1508 |        918    .0261438    .0052664         .016821    .0386513
_Ithink~1509 |        918    .0021786    .0015389         .000264    .0078477
-------------+---------------------------------------------------------------
_Ithink~1510 |        918    .0076253    .0028711        .0030711    .0156475

. 
. recode thinking (0=1) (else=0), gen(think_indep)
(21840 differences between thinking and think_indep)

. recode thinking (1501=1) (else=0), gen(think_pt)
(20867 differences between thinking and think_pt)

. recode thinking (1502=1) (else=0), gen(think_pmdb)
(20867 differences between thinking and think_pmdb)

. recode thinking (1503=1) (else=0), gen(think_psdb)
(20867 differences between thinking and think_psdb)

. recode thinking (1504=1) (else=0), gen(think_psb)
(20867 differences between thinking and think_psb)

. recode thinking (1505=1) (else=0), gen(think_dem)
(20867 differences between thinking and think_dem)

. recode thinking (1506=1) (else=0), gen(think_ptb)
(20867 differences between thinking and think_ptb)

. recode thinking (1507=1) (else=0), gen(think_psol)
(20867 differences between thinking and think_psol)

. recode thinking (1508=1) (else=0), gen(think_pv)
(20867 differences between thinking and think_pv)

. recode thinking (1509=1) (else=0), gen(think_pdt)
(20867 differences between thinking and think_pdt)

. 
. replace think_indep=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_pt=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_pmdb=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_psdb=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_psb=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_dem=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_ptb=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_psol=. if thinking==.
(20,578 real changes made, 20,578 to missing)

. replace think_pv=. if thinking==. 
(20,578 real changes made, 20,578 to missing)

. replace think_pdt=. if thinking==.
(20,578 real changes made, 20,578 to missing)

. 
. recode think_indep (1=0) (0=1), gen(think_partisan)
(1262 differences between think_indep and think_partisan)

. 
. *test-retest
. gen l5pref_pt=l5.pref_pt
(20,831 missing values generated)

. gen l5pref_psdb=l5.pref_psdb
(20,831 missing values generated)

. polychoric l5pref_pt pref_pt if wave==6

Variables :  l5pref_pt pref_pt
Type :       polychoric
Rho        = .64465676
S.e.       = .06376146
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. polychoric l5pref_psdb pref_psdb if wave==6

Variables :  l5pref_psdb pref_psdb
Type :       polychoric
Rho        = .75513635
S.e.       = .09691173
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. gen l5iden_pt=l5.iden_pt
(20,859 missing values generated)

. gen l5iden_psdb=l5.iden_psdb
(20,859 missing values generated)

. polychoric l5iden_pt iden_pt if wave==6

Variables :  l5iden_pt iden_pt
Type :       polychoric
Rho        = .56253197
S.e.       = .07698168
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. polychoric l5iden_psdb iden_psdb if wave==6

Variables :  l5iden_psdb iden_psdb
Type :       polychoric
Rho        = .42481922
S.e.       = .20047111
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. gen l5think_pt=l5.think_pt
(20,922 missing values generated)

. gen l5think_psdb=l5.think_psdb
(20,922 missing values generated)

. polychoric l5think_pt think_pt if wave==6

Variables :  l5think_pt think_pt
Type :       polychoric
Rho        = .43712745
S.e.       = .10620668
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. polychoric l5think_psdb think_psdb if wave==6

Variables :  l5think_psdb think_psdb
Type :       polychoric
Rho        = .58346169
S.e.       = .16643629
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. gen l5symp_pt=l5.symp_pt
(18,203 missing values generated)

. gen l5symp_psdb=l5.symp_psdb
(18,203 missing values generated)

. polychoric l5symp_pt symp_pt if wave==6

Variables :  l5symp_pt symp_pt
Type :       polychoric
Rho        = .6236332
S.e.       = .04323824
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. polychoric l5symp_psdb symp_psdb if wave==6

Variables :  l5symp_psdb symp_psdb
Type :       polychoric
Rho        = .71251208
S.e.       = .06946865
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. *II. OTHER MEASURES RELATED TO PARTY ID
.         *1: Party dislike 
. recode vb15 (1501=0) (else=1), gen(dislike_pt)
(21840 differences between vb15 and dislike_pt)

. replace dislike_pt=. if wave~=1
(18,720 real changes made, 18,720 to missing)

. 
. recode vb15 (1503=0) (else=1), gen(dislike_psdb)
(21840 differences between vb15 and dislike_psdb)

. replace dislike_psdb=. if wave~=1
(18,720 real changes made, 18,720 to missing)

. 
. recode vb15 (1502=0) (else=1), gen(dislike_pmdb)
(21840 differences between vb15 and dislike_pmdb)

. replace dislike_pmdb=. if wave~=1
(18,720 real changes made, 18,720 to missing)

. 
.         *2: Presidential vote choice 
. recode vb1 (1=1) (.=.) (.a=.) (.b=.) (.c=.) (else=0) , gen(vote_pt)
(4434 differences between vb1 and vote_pt)

. recode vb3 (1=1) (.=.) (.a=.) (.b=.) (.c=.) (else=0) , gen(temp)
(1110 differences between vb3 and temp)

. replace vote_pt=temp if wave==7
(972 real changes made)

. drop temp

. 
. recode vb1 (2=1) (.=.) (.a=.) (.b=.) (.c=.) (else=0) , gen(vote_psdb)
(6804 differences between vb1 and vote_psdb)

. recode vb3 (2=1) (.=.) (.a=.) (.b=.) (.c=.) (else=0) , gen(temp)
(2208 differences between vb3 and temp)

. replace vote_psdb=temp if wave==7
(972 real changes made)

. 
. 
. *Criterion-Related Validity graphs
. recode m1 1=5 2=4 3=3 4=2 5=1
(m1: 3889 changes made)

. recode vb1 1=1 2=0 .=. else=., gen(votebin)
(4434 differences between vb1 and votebin)

. 
. xtset idtelefone wave
       panel variable:  idtelefone (strongly balanced)
        time variable:  wave, 1 to 7
                delta:  1 unit

. gen l5vb1=l5.vb1
(18,915 missing values generated)

. 
. gen pref_symp=.
(21,840 missing values generated)

. replace pref_symp=-2 if symp_pt==1 
(1,169 real changes made)

. replace pref_symp=-1 if pref_pt==1 & symp==0
(150 real changes made)

. replace pref_symp=0 if preference==0 & symp==0
(819 real changes made)

. replace pref_symp=1 if pref_psdb==1 & symp==0 
(29 real changes made)

. replace pref_symp=2 if symp_psdb==1
(242 real changes made)

. replace pref_symp=. if wave>1 & wave<6
(496 real changes made, 496 to missing)

. replace pref_symp=. if wave==7
(235 real changes made, 235 to missing)

. replace pref_symp=. if pref_symp==2 & preference~=1503 & preference~=0
(79 real changes made, 79 to missing)

. replace pref_symp=. if pref_symp==-2 & preference~=1501 & preference~=0
(361 real changes made, 361 to missing)

. 
. gen dilma_minus_aecio=ft3-ft2
(19,926 missing values generated)

. bysort pref_symp: summ votebin 

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |        157    .9808917    .1373437          0          1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |         94    .9361702    .2457602          0          1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 0

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |        263    .6007605    .4906758          0          1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |          8         .25      .46291          0          1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |         29    .1724138    .3844259          0          1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = .

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |      2,695    .7283859    .4448745          0          1


. bysort pref_symp: summ dilma_minus_aecio

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |         86    3.197674    3.346006         -5         10

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |         85    2.494118    3.856422         -8         10

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 0

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |        313   -.1533546    2.823118         -8         10

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |         21   -2.619048    3.891994        -10          6

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |         26   -2.192308    3.249852         -8          4

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = .

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
dilma_minu~o |      1,383      .45047    3.562788        -10         10


. bysort pref_symp: summ m1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |        195    3.753846    .8191793          1          5

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = -1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |        149    3.308725    1.032573          1          5

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 0

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |        804    2.748756    1.018661          1          5

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |         29    2.172414    1.037475          1          4

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = 2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |         42    2.285714    1.110608          1          5

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> pref_symp = .

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
          m1 |      5,497    3.001455     1.08617          1          5


. 
. xtset idtelefone wave
       panel variable:  idtelefone (strongly balanced)
        time variable:  wave, 1 to 7
                delta:  1 unit

. gen l5pref_symp=l5.pref_symp
(20,975 missing values generated)

. capture label drop party

. label define party 2 "Agreed PSDBistas" 1 "Cont'ed PSDBistas" 0 "Agreed nonpart's" -1 "Cont'ed PTistas" -2 "Agreed PTistas"

. label values pref_symp party

. label values l5pref_symp party

. 
. gen votefuture=votebin
(18,594 missing values generated)

. replace votefuture=0 if l5pref_symp>0
(20,208 real changes made)

. 
. recode pref_symp -1=1 1=1 0=2 .=. else=0, gen(fns)
(1209 differences between pref_symp and fns)

. recode l5pref_symp -1=1 1=1 0=2 .=. else=0, gen(fns2)
(836 differences between l5pref_symp and fns2)

. separate votebin, by(fns)

              storage   display    value
variable name   type    format     label      variable label
-------------------------------------------------------------------------------------------------------------------------------------------------------------
votebin0        byte    %9.0g                 votebin, fns == 0
votebin1        byte    %9.0g                 votebin, fns == 1
votebin2        byte    %9.0g                 votebin, fns == 2

. separate votefuture, by(fns2)

              storage   display    value
variable name   type    format     label      variable label
-------------------------------------------------------------------------------------------------------------------------------------------------------------
votefuture0     byte    %9.0g                 votefuture, fns2 == 0
votefuture1     byte    %9.0g                 votefuture, fns2 == 1
votefuture2     byte    %9.0g                 votefuture, fns2 == 2

. separate m1, by(fns)

              storage   display    value
variable name   type    format     label      variable label
-------------------------------------------------------------------------------------------------------------------------------------------------------------
m10             byte    %8.0g      M1         m1, fns == 0
m11             byte    %8.0g      M1         m1, fns == 1
m12             byte    %8.0g      M1         m1, fns == 2

. separate dilma_minus_aecio, by(fns)

              storage   display    value
variable name   type    format     label      variable label
-------------------------------------------------------------------------------------------------------------------------------------------------------------
dilma_minus_a~0 byte    %9.0g                 dilma_minus_aecio, fns == 0
dilma_minus_a~1 byte    %9.0g                 dilma_minus_aecio, fns == 1
dilma_minus_a~2 byte    %9.0g                 dilma_minus_aecio, fns == 2

. 
.         *bar graphs
. graph bar (mean) votebin?, over(pref_symp, label(angle(0) labsize(small) alternate tick labgap(0)))  graphregion(color(white)) plotregion(lcolor(black)) //
> /
> blabel(bar, format(%9.2f) size(small)) bar(1, color(black)) bar(2, color(gs12)) bar(3, color(gs16) lcolor(black)) title("A1: Rousseff current 2-party vote"
> ) ytitle("") ylab(, labsize(small))  legend(off) nofill

. graph save "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr1.gph", replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\gr1.gph saved)

. 
. graph bar (mean) votefuture? if wave==6, over(l5pref_symp, label(angle(0) labsize(small) alternate tick labgap(0)))  graphregion(color(white)) plotregion(l
> color(black)) ///
> blabel(bar, format(%9.2f) size(small)) bar(3, color(gs16) lcolor(black)) bar(2, color(gs12)) bar(1, color(black)) title("A2: Rousseff future 2-party vote")
>  ytitle("") ylab(, labsize(small)) text(.05 81 "Not enough observations", size(small) box fcolor(gs12) lcolor(white)) legend(off) nofill

. graph save "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr2.gph", replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\gr2.gph saved)

. 
. graph bar (mean) m1?, over(pref_symp, label(angle(0) labsize(small) alternate tick labgap(0))) graphregion(color(white)) plotregion(lcolor(black))  ///
> blabel(bar, format(%9.2f) size(small)) bar(1, color(black)) bar(2, color(gs12)) bar(3, color(gs16) lcolor(black))   title("A3: Rousseff job approval") ytit
> le("") exclude0 yscale(range(2(.5)4)) ytick(2 2.5 3 3.5 4) ylab(2 2.5 3 3.5 4)  ylab(, labsize(small))  legend(off) nofill

. graph save "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr3.gph", replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\gr3.gph saved)

. 
. graph bar (mean) dilma_minus_aecio?, over(pref_symp, label(angle(0) labsize(small) alternate tick labgap(0))) graphregion(color(white)) plotregion(lcolor(b
> lack))  ///
> blabel(bar, format(%9.2f) size(small)) bar(1, color(black)) bar(2, color(gs12)) bar(3, color(gs16) lcolor(black))   title("A4: Rousseff - Neves feeling the
> rms") ytitle("") ytick(-5 -2.5 0 2.5 5) ylab(-5 -2.5 0 2.5 5)  ylab(, labsize(small))  legend(off) nofill

. graph save "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr4.gph", replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\gr4.gph saved)

. 
. *local stick f
. graph combine ///
> "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr1.gph" ///
> "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr2.gph" ///
> "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr3.gph" ///
> "`stick'\My Documents\Research\Paper PID Experiment\Brazil\gr4.gph", ///
> row(2) col(2) graphregion(color(white))

. graph export "`stick'\My Documents\Research\Paper PID Experiment\LaTeX\JOPBrazilCriterion.pdf", as(pdf) replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\LaTeX\JOPBrazilCriterion.pdf wri
> tten in PDF format)

. 
. *III. Descriptives
. summ pref_pt iden_pt think_pt symp_pt dislike_pt vote_pt ft6 pref_psdb iden_psdb think_psdb symp_psdb dislike_psdb vote_psdb ft7 if wave==1 | wave==6

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     pref_pt |      1,418     .248237    .4321424          0          1
     iden_pt |      1,409    .1760114     .380965          0          1
    think_pt |      1,262     .118859    .3237508          0          1
     symp_pt |      4,245    .1316843    .3381871          0          1
  dislike_pt |      3,120    .8541667    .3529956          0          1
-------------+---------------------------------------------------------
     vote_pt |      3,653    .3939228     .488685          0          1
         ft6 |      2,828    4.946252    3.443815          0         10
   pref_psdb |      1,418    .0521862    .2224808          0          1
   iden_psdb |      1,409    .0390348    .1937465          0          1
  think_psdb |      1,262    .0285261    .1665364          0          1
-------------+---------------------------------------------------------
   symp_psdb |      4,245    .0285041    .1664276          0          1
dislike_psdb |      3,120    .9557692    .2056404          0          1
   vote_psdb |      3,653    .1554886    .3624195          0          1
         ft7 |      2,416     4.15894    2.819994          0         10

. summ votebin m1 dilma_minus_aecio if pref_symp~=.

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |        551     .738657    .4397657          0          1
          m1 |      1,219    2.948318    1.078124          1          5
dilma_minu~o |        531    .6158192     3.56546        -10         10

. summ votebin if l5pref_symp~=. & wave==6

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     votebin |        200        .695    .4615628          0          1

. 
. tab iden_pt vote_pt if wave==1, row

+----------------+
| Key            |
|----------------|
|   frequency    |
| row percentage |
+----------------+

 RECODE of |
  identity |
(RECODE of | RECODE of vb1 (v1. em
exp_2 (exp |    outubro havera
2. voce se |     eleicoes para
 considera |   presidente. se a
pedetista, |     eleicao foss
   pemedeb |         0          1 |     Total
-----------+----------------------+----------
         0 |       448        184 |       632 
           |     70.89      29.11 |    100.00 
-----------+----------------------+----------
         1 |        33        124 |       157 
           |     21.02      78.98 |    100.00 
-----------+----------------------+----------
     Total |       481        308 |       789 
           |     60.96      39.04 |    100.00 


. tab think_pt vote_pt if wave==1, row

+----------------+
| Key            |
|----------------|
|   frequency    |
| row percentage |
+----------------+

           | RECODE of vb1 (v1. em
           |    outubro havera
           |     eleicoes para
           |   presidente. se a
 RECODE of |     eleicao foss
  thinking |         0          1 |     Total
-----------+----------------------+----------
         0 |       462        191 |       653 
           |     70.75      29.25 |    100.00 
-----------+----------------------+----------
         1 |        11         75 |        86 
           |     12.79      87.21 |    100.00 
-----------+----------------------+----------
     Total |       473        266 |       739 
           |     64.01      35.99 |    100.00 


. tab pref_pt vote_pt if wave==1, row 

+----------------+
| Key            |
|----------------|
|   frequency    |
| row percentage |
+----------------+

 RECODE of |
preference |
(RECODE of |
exp_1 (exp | RECODE of vb1 (v1. em
 1. qual e |    outubro havera
     o seu |     eleicoes para
   partido |   presidente. se a
  politico |     eleicao foss
      de p |         0          1 |     Total
-----------+----------------------+----------
         0 |       445        105 |       550 
           |     80.91      19.09 |    100.00 
-----------+----------------------+----------
         1 |        66        184 |       250 
           |     26.40      73.60 |    100.00 
-----------+----------------------+----------
     Total |       511        289 |       800 
           |     63.88      36.13 |    100.00 


. tab symp_pt vote_pt if wave==1, row 

+----------------+
| Key            |
|----------------|
|   frequency    |
| row percentage |
+----------------+

 RECODE of |
      vb11 |
(vb11. com |
      qual | RECODE of vb1 (v1. em
   partido |    outubro havera
      o(a) |     eleicoes para
   sr.(a). |   presidente. se a
simpatiza? |     eleicao foss
         ) |         0          1 |     Total
-----------+----------------------+----------
         0 |     1,463        638 |     2,101 
           |     69.63      30.37 |    100.00 
-----------+----------------------+----------
         1 |        59        262 |       321 
           |     18.38      81.62 |    100.00 
-----------+----------------------+----------
     Total |     1,522        900 |     2,422 
           |     62.84      37.16 |    100.00 


. 
.  
. *IV. Randomization check
. *income, but not pretreatment
. *recode p15 (1=724) (2=1810) (3=2896) (4=4711) (5=6516) (6=10860) (7=18100)
. *replace p15=ln(p15)
. recode p1 2=1 1=0
(p1: 3120 changes made)

. recode urb_rur 2=0
(urb_rur: 821 changes made)

. bysort vb10exp: ci means symp_indep symp_pt p1 idade pol1 urb_rur if wave==1

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> vb10exp = 1

    Variable |        Obs        Mean    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
  symp_indep |      1,029    .7793975    .0129327          .75402    .8047749
     symp_pt |      1,029    .1311953    .0105299        .1105328    .1518579
          p1 |      1,058          .5    .0153792        .4698229    .5301771
       idade |      1,058    40.15028    .5180663        39.13373    41.16684
        pol1 |      1,055    2.989573    .0303889        2.929944    3.049203
-------------+---------------------------------------------------------------
     urb_rur |      1,058    .8742911     .010197        .8542824    .8942998

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> vb10exp = 2

    Variable |        Obs        Mean    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
  symp_indep |      1,032    .7945736    .0125825        .7698835    .8192638
     symp_pt |      1,032    .1133721     .009874        .0939966    .1327476
          p1 |      1,060    .5056604    .0153636        .4755137     .535807
       idade |      1,060    40.10283     .515992        39.09035    41.11531
        pol1 |      1,057    3.005676    .0298712        2.947063     3.06429
-------------+---------------------------------------------------------------
     urb_rur |      1,060    .9509434    .0066371          .93792    .9639667

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> vb10exp = 3

    Variable |        Obs        Mean    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
  symp_indep |        981    .8103976    .0125215        .7858254    .8349697
     symp_pt |        981    .1090724    .0099578        .0895312    .1286135
          p1 |      1,002      .50499    .0158027        .4739798    .5360002
       idade |      1,002     39.8982    .5223433        38.87319    40.92322
        pol1 |      1,000       3.038    .0292132        2.980674    3.095326
-------------+---------------------------------------------------------------
     urb_rur |      1,002    .8512974    .0112456        .8292297    .8733651

-------------------------------------------------------------------------------------------------------------------------------------------------------------
-> vb10exp = .

    Variable |        Obs        Mean    Std. Err.       [95% Conf. Interval]
-------------+---------------------------------------------------------------
  symp_indep |          0           .           .               .           .
     symp_pt |          0           .           .               .           .
          p1 |          0           .           .               .           .
       idade |          0           .           .               .           .
        pol1 |          0           .           .               .           .
-------------+---------------------------------------------------------------
     urb_rur |          0           .           .               .           .

. 
. *V: Factor Analysis check for unidimensionality
. polychoricpca pref_pt symp_pt vote_pt ft6

Polychoric correlation matrix

           pref_pt    symp_pt    vote_pt        ft6
pref_pt          1
symp_pt  .89030238          1
vote_pt  .79009313  .68982375          1
    ft6  .74745321  .65773646  .67499872          1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    3.230814   |    0.807703            |   0.807703
  2 |    0.368201   |    0.092050            |   0.899754
  3 |    0.312774   |    0.078194            |   0.977947
  4 |    0.088211   |    0.022053            |   1.000000

. polychoricpca iden_pt symp_pt vote_pt ft6

Polychoric correlation matrix

           iden_pt    symp_pt    vote_pt        ft6
iden_pt          1
symp_pt  .82789673          1
vote_pt  .67667856  .60970741          1
    ft6  .71454002  .71698506  .61977547          1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    3.087972   |    0.771993            |   0.771993
  2 |    0.427152   |    0.106788            |   0.878781
  3 |    0.320003   |    0.080001            |   0.958782
  4 |    0.164873   |    0.041218            |   1.000000

. polychoricpca think_pt symp_pt vote_pt ft6

Polychoric correlation matrix

           think_pt    symp_pt    vote_pt        ft6
think_pt          1
 symp_pt   .7744154          1
 vote_pt  .73171622  .70596549          1
     ft6  .74243845  .66153529  .67324923          1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    3.146079   |    0.786520            |   0.786520
  2 |    0.346200   |    0.086550            |   0.873070
  3 |    0.302282   |    0.075570            |   0.948640
  4 |    0.205439   |    0.051360            |   1.000000

. 
. polychoricpca pref_psdb symp_psdb vote_psdb ft7

Polychoric correlation matrix

           pref_psdb  symp_psdb  vote_psdb        ft7
pref_psdb          1
symp_psdb  .96436592          1
vote_psdb  .58408038  .62229871          1
      ft7  .62274337  .66401416  .29741848          1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    2.921955   |    0.730489            |   0.730489
  2 |    0.703773   |    0.175943            |   0.906432
  3 |    0.342085   |    0.085521            |   0.991953
  4 |    0.032186   |    0.008047            |   1.000000

. polychoricpca iden_psdb symp_psdb vote_psdb ft7

Polychoric correlation matrix

           iden_psdb  symp_psdb  vote_psdb        ft7
iden_psdb          1
symp_psdb  .90008076          1
vote_psdb  .52868729  .59519667          1
      ft7  .54888214  .68976952  .34863813          1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    2.839855   |    0.709964            |   0.709964
  2 |    0.656708   |    0.164177            |   0.874141
  3 |    0.428162   |    0.107041            |   0.981181
  4 |    0.075274   |    0.018819            |   1.000000

. polychoricpca think_psdb symp_psdb vote_psdb ft7

Polychoric correlation matrix

            think_psdb   symp_psdb   vote_psdb         ft7
think_psdb           1
 symp_psdb   .81685592           1
 vote_psdb   .47286304   .51852087           1
       ft7   .53315016   .66591531   .35888353           1

Principal component analysis

 k  |  Eigenvalues  |  Proportion explained  |  Cum. explained
----+---------------+------------------------+------------------
  1 |    2.711243   |    0.677811            |   0.677811
  2 |    0.660966   |    0.165242            |   0.843052
  3 |    0.468228   |    0.117057            |   0.960109
  4 |    0.159563   |    0.039891            |   1.000000

. 
. *VI. IRT models
. *Note: I often fed starting values to these IRT models but it never seemed to make a difference. What does make a difference is number of integration point
> s and the integration method.
. *More intpoints are more accurate (Ayala 2009, p. 75 and 104), and ghermite is not trustworthy (see "semintro12" Stata manual, p. 6). Models using ghermite
>  would often converge but were very unstable/sensitive and often nonsensical SEs. 
. *mcaghermite produced stable results repeatedly, although slow. mvaghermite also stable and trustworthy, but slightly less likely to converge
. *Small parties are more challenging to estimate because of skew (Ayala 2009, 104), so I did these with 60 intpoints
. 
.         *1: PT Results
. irt grm pref_pt iden_pt think_pt symp_pt vote_pt ft6 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmethod(mcaghermite) 

Fitting fixed-effects model:

Iteration 0:   log likelihood = -12422.421  
Iteration 1:   log likelihood = -12422.421  

Fitting full model:

Iteration 0:   log pseudolikelihood = -11659.226  
Iteration 1:   log pseudolikelihood = -10913.224  
Iteration 2:   log pseudolikelihood = -10866.138  
Iteration 3:   log pseudolikelihood = -10838.203  
Iteration 4:   log pseudolikelihood = -10836.125  
Iteration 5:   log pseudolikelihood = -10835.961  
Iteration 6:   log pseudolikelihood =  -10835.97  
Iteration 7:   log pseudolikelihood = -10835.968  
Iteration 8:   log pseudolikelihood = -10835.968  

Graded response model                           Number of obs     =      4,325
Log pseudolikelihood = -10835.968
                         (Std. Err. adjusted for 3,119 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_pt      |
     Discrim |   8.426108   1.550727     5.43   0.000     5.386739    11.46548
        Diff |
         =1  |   .7293083   .0340742    21.40   0.000     .6625241    .7960926
-------------+----------------------------------------------------------------
iden_pt      |
     Discrim |   5.282047   1.079013     4.90   0.000      3.16722    7.396873
        Diff |
         =1  |   1.004126   .0416328    24.12   0.000      .922527    1.085725
-------------+----------------------------------------------------------------
think_pt     |
     Discrim |   4.742486    .821964     5.77   0.000     3.131466    6.353506
        Diff |
         =1  |   1.266129   .0523903    24.17   0.000     1.163445    1.368812
-------------+----------------------------------------------------------------
symp_pt      |
     Discrim |   4.104394   .3066927    13.38   0.000     3.503288    4.705501
        Diff |
         =1  |   1.227347   .0340247    36.07   0.000     1.160659    1.294034
-------------+----------------------------------------------------------------
vote_pt      |
     Discrim |   1.930926    .103649    18.63   0.000     1.727778    2.134075
        Diff |
         =1  |   .3970237   .0323799    12.26   0.000     .3335603    .4604871
-------------+----------------------------------------------------------------
ft6          |
     Discrim |   2.233214   .0969018    23.05   0.000      2.04329    2.423139
        Diff |
       >= 1  |  -1.086701   .0404161   -26.89   0.000    -1.165915   -1.007487
       >= 2  |  -.9258942   .0379385   -24.41   0.000    -1.000252   -.8515361
       >= 3  |  -.6992054   .0348161   -20.08   0.000    -.7674438   -.6309671
       >= 4  |  -.4825819   .0325779   -14.81   0.000    -.5464335   -.4187303
       >= 5  |  -.3026303   .0311602    -9.71   0.000    -.3637031   -.2415575
       >= 6  |   .1855547    .029623     6.26   0.000     .1274947    .2436146
       >= 7  |   .3974988   .0301693    13.18   0.000      .338368    .4566296
       >= 8  |   .7103959   .0326911    21.73   0.000     .6463225    .7744692
       >= 9  |   1.206417   .0409062    29.49   0.000     1.126242    1.286591
        =10  |   1.486331   .0483261    30.76   0.000     1.391614    1.581049
------------------------------------------------------------------------------

. test _b[pref_pt:Theta]=_b[think_pt:Theta]

 ( 1)  [pref_pt]Theta - [think_pt]Theta = 0

           chi2(  1) =    4.44
         Prob > chi2 =    0.0350

. test _b[pref_pt:Theta]=_b[symp_pt:Theta]

 ( 1)  [pref_pt]Theta - [symp_pt]Theta = 0

           chi2(  1) =    7.49
         Prob > chi2 =    0.0062

. test _b[pref_pt:Theta]=_b[iden_pt:Theta]

 ( 1)  [pref_pt]Theta - [iden_pt]Theta = 0

           chi2(  1) =    2.81
         Prob > chi2 =    0.0938

. predict thetapt, latent
(option ebmeans assumed)
(using 30 quadrature points)

. summ thetapt

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     thetapt |     21,840    .0222707    .4707131  -1.198563   1.937335

. 
. irtgraph icc (ft6 , lcolor(gs14) lwidth(vthin) lpattern(dash))  (vote_pt, lcolor(black) lwidth(medium) lpattern(dash)) (pref_pt iden_pt symp_pt think_pt, l
> color(black) lwidth(medthick)), ///
> addplot((scatteri .5 0.33 "Pres. Vote", msymbol(i) mlabpos(0) mlabangle(52) mlabcolor(black) mlabsize(vsmall)) ///
> (scatteri .5 .67 "Prefer({it:P{superscript:0}N{superscript:0}})", msymbol(i) mlabpos(0) mlabangle(82) mlabcolor(black) mlabsize(small)) ///
> (scatteri .5 .94 "Consider({it:P{superscript:1}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(74) mlabcolor(black) mlabsize(small)) ///
> (scatteri .4 1.06 "Sympathy({it:P{superscript:0}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(69) mlabcolor(black) mlabsize(small)) /// 
> (scatteri .4 1.27 "Way of Thinking({it:P{superscript:0}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(71) mlabcolor(black) mlabsize(small))) ///
> bcc range(-1.20 1.94) legend(off) graphregion(color(white)) plotregion(lstyle(yxline) lcolor(black)) b1title("{it: θ{subscript:PT}}: Strength of PT Partisa
> nship") xtitle("") ytitle("{it:Pr}({it:x{subscript: j•PT}}|{it: θ{subscript:PT}})")  title("") xtick(-2(.5)2) xlabel(-2(1)2) 

. graph export "`stick'\My Documents\Research\Paper PID Experiment\LaTeX\JOPPTIRT.pdf", as(pdf) replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\LaTeX\JOPPTIRT.pdf written in PD
> F format)

. 
. polychoric thetapt pref_pt iden_pt think_pt symp_pt vote_pt ft6 if wave==1 | wave==6, pw verbose

Variables :  thetapt pref_pt
Type :       polyserial
Rho        = .93383082
S.e.       = .0120561

Variables :  thetapt iden_pt
Type :       polyserial
Rho        = .76109127
S.e.       = .01785293

Variables :  thetapt think_pt
Type :       polyserial
Rho        = .64631043
S.e.       = .01790274

Variables :  thetapt symp_pt
Type :       polyserial
Rho        = .61836566
S.e.       = .00938197

Variables :  thetapt vote_pt
Type :       polyserial
Rho        = .6250139
S.e.       = .01327066

Variables :  thetapt ft6
Type :       Pearson
Rho        = .89003733
S.e.       = .00857574
no observations

Variables :  pref_pt iden_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt symp_pt
Type :       polychoric
Rho        = .9187894
S.e.       = .01423103
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt vote_pt
Type :       polychoric
Rho        = .71012471
S.e.       = .02974057
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt ft6
Type :       polyserial
Rho        = .75496226
S.e.       = .02056818
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt symp_pt
Type :       polychoric
Rho        = .87293879
S.e.       = .02002008
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt vote_pt
Type :       polychoric
Rho        = .64405588
S.e.       = .03632682
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt ft6
Type :       polyserial
Rho        = .70861935
S.e.       = .03472193

Variables :  think_pt symp_pt
Type :       polychoric
Rho        = .86335961
S.e.       = .02359722
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt vote_pt
Type :       polychoric
Rho        = .67232313
S.e.       = .04102964
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt ft6
Type :       polyserial
Rho        = .71945907
S.e.       = .03809554

Variables :  symp_pt vote_pt
Type :       polychoric
Rho        = .65712393
S.e.       = .02213691
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  symp_pt ft6
Type :       polyserial
Rho        = .6735074
S.e.       = .02193122

Variables :  vote_pt ft6
Type :       polyserial
Rho        = .65456931
S.e.       = .01560988

.         
.         *2: PSDB Results
. set more off

. *to speed things, I did feed starting values here.  They came from this command, which took several hours to converge:
. *irt grm pref_psdb iden_psdb think_psdb symp_psdb dislike_psdb vote_psdb ft7 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmetho
> d(mcaghermite) 
. mat start = [7.7449082, 3.6025677, 3.4145009, 8.536266, 1.1026588, 1.0339763, 1.338271, 13.02818, 7.0484067, 7.3877823, 16.618127, -3.6155326, 2.0287019, -
> 1.9292317, -1.656367, -1.1379195, -.62497109, -.09281607, 1.0717178, 1.7532436, 2.5152163, 3.6803523, 4.2057435, 1]

. mat colnames start = pref_psdb:Theta iden_psdb:Theta think_psdb:Theta symp_psdb:Theta dislike_psdb:Theta vote_psdb:Theta ft7:Theta pref_psdb_cut1:_cons ide
> n_psdb_cut1:_cons think_psdb_cut1:_cons symp_psdb_cut1:_cons dislike_psdb_cut1:_cons vote_psdb_cut1:_cons ft7_cut1:_cons ft7_cut2:_cons ft7_cut3:_cons ft7_
> cut4:_cons ft7_cut5:_cons ft7_cut6:_cons ft7_cut7:_cons ft7_cut8:_cons ft7_cut9:_cons ft7_cut10:_cons var(Theta):_cons

. 
. display "$S_TIME  $S_DATE"
14:06:16   6 Mar 2018

. irt grm pref_psdb iden_psdb think_psdb symp_psdb vote_psdb ft7 if wave==1 | wave==6, intpoints(60) difficult vce(cluster idtelefone) intmethod(mcaghermite)
>  from(start, skip)

Iteration 0:   log pseudolikelihood = -7722.3763  
Iteration 1:   log pseudolikelihood = -7722.1548  
Iteration 2:   log pseudolikelihood = -7722.1411  
Iteration 3:   log pseudolikelihood = -7722.1408  

Graded response model                           Number of obs     =      4,325
Log pseudolikelihood = -7722.1408
                         (Std. Err. adjusted for 3,119 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_psdb    |
     Discrim |   6.969716   1.834195     3.80   0.000     3.374759    10.56467
        Diff |
         =1  |   1.695684   .0489386    34.65   0.000     1.599766    1.791602
-------------+----------------------------------------------------------------
iden_psdb    |
     Discrim |   3.738033   .7227442     5.17   0.000     2.321481    5.154586
        Diff |
         =1  |   1.936836    .097962    19.77   0.000     1.744834    2.128838
-------------+----------------------------------------------------------------
think_psdb   |
     Discrim |   3.373295   .6537058     5.16   0.000     2.092055    4.654534
        Diff |
         =1  |    2.16365   .1330118    16.27   0.000     1.902951    2.424348
-------------+----------------------------------------------------------------
symp_psdb    |
     Discrim |   10.01014   3.508205     2.85   0.004     3.134187     16.8861
        Diff |
         =1  |   1.926694   .0700932    27.49   0.000     1.789314    2.064074
-------------+----------------------------------------------------------------
vote_psdb    |
     Discrim |   1.050479   .0940618    11.17   0.000     .8661212    1.234837
        Diff |
         =1  |   1.940698   .1369605    14.17   0.000      1.67226    2.209136
-------------+----------------------------------------------------------------
ft7          |
     Discrim |   1.317125   .0989573    13.31   0.000     1.123172    1.511078
        Diff |
       >= 1  |  -1.441507   .0842393   -17.11   0.000    -1.606613   -1.276401
       >= 2  |  -1.236659   .0753748   -16.41   0.000    -1.384391   -1.088927
       >= 3  |  -.8476585   .0600692   -14.11   0.000    -.9653919   -.7299251
       >= 4  |  -.4624662   .0481572    -9.60   0.000    -.5568526   -.3680798
       >= 5  |  -.0624777   .0415647    -1.50   0.133     -.143943    .0189877
       >= 6  |   .8138863   .0554888    14.67   0.000     .7051302    .9226424
       >= 7  |   1.327365   .0755192    17.58   0.000      1.17935     1.47538
       >= 8  |   1.901312   .1040442    18.27   0.000      1.69739    2.105235
       >= 9  |   2.780345    .160458    17.33   0.000     2.465853    3.094837
        =10  |   3.177681   .1903897    16.69   0.000     2.804524    3.550838
------------------------------------------------------------------------------

. display "$S_TIME  $S_DATE"
15:15:06   6 Mar 2018

. 
. test _b[pref_psdb:Theta]=_b[think_psdb:Theta]

 ( 1)  [pref_psdb]Theta - [think_psdb]Theta = 0

           chi2(  1) =    3.43
         Prob > chi2 =    0.0641

. test _b[pref_psdb:Theta]=_b[symp_psdb:Theta]

 ( 1)  [pref_psdb]Theta - [symp_psdb]Theta = 0

           chi2(  1) =    0.48
         Prob > chi2 =    0.4867

. test _b[pref_psdb:Theta]=_b[iden_psdb:Theta]

 ( 1)  [pref_psdb]Theta - [iden_psdb]Theta = 0

           chi2(  1) =    2.62
         Prob > chi2 =    0.1053

. predict thetapsdb, latent 
(option ebmeans assumed)
(using 60 quadrature points)

. summ thetapsdb

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
   thetapsdb |     21,840    .0062158    .3360292  -.9148247   2.630779

. 
. irtgraph icc (ft7, lcolor(gs14) lwidth(vthin) lpattern(dash))  (vote_psdb, lcolor(black) lwidth(medium) lpattern(dash)) (pref_psdb iden_psdb symp_psdb thin
> k_psdb, lcolor(black) lwidth(medthick)), ///
> addplot((scatteri .375 1.37 "Pres. Vote", msymbol(i) mlabpos(0) mlabangle(37) mlabcolor(black) mlabsize(vsmall)) ///
> (scatteri .55 1.65 "Prefer({it:P{superscript:0}N{superscript:0}})", msymbol(i) mlabpos(0) mlabangle(82) mlabcolor(black) mlabsize(small)) ///
> (scatteri .88 2.39 "Consider({it:P{superscript:1}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(55) mlabcolor(black) mlabsize(small)) ///
> (scatteri .62 1.9 "Sympathy({it:P{superscript:0}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(83) mlabcolor(black) mlabsize(small)) /// 
> (scatteri .4 2.14 "Way of Thinking({it:P{superscript:0}N{superscript:1}})", msymbol(i) mlabpos(0) mlabangle(67) mlabcolor(black) mlabsize(small))) ///
> bcc range(-.92 2.64) legend(off) graphregion(color(white)) plotregion(lstyle(yxline) lcolor(black)) b1title("{it: θ{subscript:PSDB}}: Strength of PSDB Part
> isanship") xtitle("") ytitle("{it:Pr}({it:x{subscript: j•PSDB}}|{it: θ{subscript:PSDB}})") title("") xtick(-2(.5)3) xlabel(-2(1)3)

. graph export "`stick'\My Documents\Research\Paper PID Experiment\LaTeX\JOPPSDBIRT.pdf", as(pdf) replace
(file C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\LaTeX\JOPPSDBIRT.pdf written in 
> PDF format)

. 
. polychoric thetapsdb pref_psdb iden_psdb think_psdb symp_psdb vote_psdb ft7 if wave==1 | wave==6, pw verbose

Variables :  thetapsdb pref_psdb
Type :       polyserial
Rho        = .49892914
S.e.       = .02262668

Variables :  thetapsdb iden_psdb
Type :       polyserial
Rho        = .42676545
S.e.       = .02577502

Variables :  thetapsdb think_psdb
Type :       polyserial
Rho        = .35706167
S.e.       = .02549694

Variables :  thetapsdb symp_psdb
Type :       polyserial
Rho        = .39140146
S.e.       = .01048681

Variables :  thetapsdb vote_psdb
Type :       polyserial
Rho        = .38457642
S.e.       = .01848584

Variables :  thetapsdb ft7
Type :       Pearson
Rho        = .84828217
S.e.       = .01077789
no observations

Variables :  pref_psdb iden_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb symp_psdb
Type :       polychoric
Rho        = .961506
S.e.       = .01442384
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb vote_psdb
Type :       polychoric
Rho        = .48988055
S.e.       = .06569833
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb ft7
Type :       polyserial
Rho        = .60048799
S.e.       = .04778163
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb symp_psdb
Type :       polychoric
Rho        = .89534667
S.e.       = .03203524
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb vote_psdb
Type :       polychoric
Rho        = .44076549
S.e.       = .07467616
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb ft7
Type :       polyserial
Rho        = .53779083
S.e.       = .06446141

Variables :  think_psdb symp_psdb
Type :       polychoric
Rho        = .84448304
S.e.       = .04660182
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_psdb vote_psdb
Type :       polychoric
Rho        = .54805669
S.e.       = .07433876
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_psdb ft7
Type :       polyserial
Rho        = .53918703
S.e.       = .08239575

Variables :  symp_psdb vote_psdb
Type :       polychoric
Rho        = .51515745
S.e.       = .04299682
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  symp_psdb ft7
Type :       polyserial
Rho        = .64907714
S.e.       = .03258366

Variables :  vote_psdb ft7
Type :       polyserial
Rho        = .32811693
S.e.       = .03125678

. */
. 
. *Item-rest correlation
.         *1: PT Results
. capture drop thetapt

. irt grm iden_pt think_pt symp_pt vote_pt ft6 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmethod(mcaghermite) 

Fitting fixed-effects model:

Iteration 0:   log likelihood = -11627.789  
Iteration 1:   log likelihood = -11627.789  

Fitting full model:

Iteration 0:   log pseudolikelihood = -11051.973  
Iteration 1:   log pseudolikelihood = -10474.801  
Iteration 2:   log pseudolikelihood = -10454.765  
Iteration 3:   log pseudolikelihood = -10433.773  
Iteration 4:   log pseudolikelihood = -10433.108  
Iteration 5:   log pseudolikelihood = -10433.092  
Iteration 6:   log pseudolikelihood = -10433.092  

Graded response model                           Number of obs     =      4,323
Log pseudolikelihood = -10433.092
                         (Std. Err. adjusted for 3,117 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
iden_pt      |
     Discrim |   5.332216    1.15078     4.63   0.000     3.076728    7.587704
        Diff |
         =1  |   .9958612   .0424658    23.45   0.000     .9126298    1.079093
-------------+----------------------------------------------------------------
think_pt     |
     Discrim |   4.751191   .8654243     5.49   0.000     3.054991    6.447392
        Diff |
         =1  |   1.260399   .0530617    23.75   0.000       1.1564    1.364398
-------------+----------------------------------------------------------------
symp_pt      |
     Discrim |   3.567347   .2858961    12.48   0.000        3.007    4.127693
        Diff |
         =1  |   1.253653   .0370765    33.81   0.000     1.180985    1.326322
-------------+----------------------------------------------------------------
vote_pt      |
     Discrim |   2.110018   .1267424    16.65   0.000     1.861607    2.358428
        Diff |
         =1  |   .3814158   .0308247    12.37   0.000     .3210005    .4418311
-------------+----------------------------------------------------------------
ft6          |
     Discrim |   2.235091   .1111386    20.11   0.000     2.017263    2.452918
        Diff |
       >= 1  |  -1.100352   .0413914   -26.58   0.000    -1.181478   -1.019226
       >= 2  |   -.941209   .0386353   -24.36   0.000    -1.016933   -.8654852
       >= 3  |  -.7167817   .0351949   -20.37   0.000    -.7857625   -.6478009
       >= 4  |  -.5017763   .0326643   -15.36   0.000    -.5657971   -.4377554
       >= 5  |  -.3230793   .0310692   -10.40   0.000    -.3839738   -.2621848
       >= 6  |    .161402   .0292154     5.52   0.000     .1041408    .2186633
       >= 7  |   .3722688     .02977    12.50   0.000     .3139207     .430617
       >= 8  |   .6854316   .0325379    21.07   0.000     .6216584    .7492047
       >= 9  |   1.183585   .0414436    28.56   0.000     1.102357    1.264813
        =10  |   1.464604   .0496246    29.51   0.000     1.367341    1.561866
------------------------------------------------------------------------------

. predict thetapt, latent
(option ebmeans assumed)
(using 30 quadrature points)

. polychoric thetapt pref_pt iden_pt think_pt symp_pt, pw verbose 

Variables :  thetapt pref_pt
Type :       polyserial
Rho        = .61770338
S.e.       = .01787143

Variables :  thetapt iden_pt
Type :       polyserial
Rho        = .76312555
S.e.       = .01736675

Variables :  thetapt think_pt
Type :       polyserial
Rho        = .64600183
S.e.       = .01750577

Variables :  thetapt symp_pt
Type :       polyserial
Rho        = .67610614
S.e.       = .00691692
no observations

Variables :  pref_pt iden_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt symp_pt
Type :       polychoric
Rho        = .9187894
S.e.       = .01423103
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt symp_pt
Type :       polychoric
Rho        = .87293879
S.e.       = .02002008
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt symp_pt
Type :       polychoric
Rho        = .86335961
S.e.       = .02359723
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. drop thetapt

. irt grm pref_pt think_pt symp_pt vote_pt ft6 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmethod(mcaghermite) 

Fitting fixed-effects model:

Iteration 0:   log likelihood = -11766.826  
Iteration 1:   log likelihood = -11766.826  

Fitting full model:

Iteration 0:   log pseudolikelihood = -11140.389  
Iteration 1:   log pseudolikelihood = -10515.771  
Iteration 2:   log pseudolikelihood = -10495.341  
Iteration 3:   log pseudolikelihood = -10462.304  
Iteration 4:   log pseudolikelihood = -10460.323  
Iteration 5:   log pseudolikelihood = -10460.176  
Iteration 6:   log pseudolikelihood = -10460.177  
Iteration 7:   log pseudolikelihood = -10460.176  

Graded response model                           Number of obs     =      4,321
Log pseudolikelihood = -10460.176
                         (Std. Err. adjusted for 3,117 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_pt      |
     Discrim |   8.423207   1.674873     5.03   0.000     5.140517     11.7059
        Diff |
         =1  |    .722889   .0338215    21.37   0.000     .6566001    .7891779
-------------+----------------------------------------------------------------
think_pt     |
     Discrim |   4.770097   .8532799     5.59   0.000     3.097699    6.442495
        Diff |
         =1  |   1.263432   .0527144    23.97   0.000     1.160114     1.36675
-------------+----------------------------------------------------------------
symp_pt      |
     Discrim |   3.776762   .3138769    12.03   0.000     3.161574    4.391949
        Diff |
         =1  |   1.242736   .0365107    34.04   0.000     1.171176    1.314295
-------------+----------------------------------------------------------------
vote_pt      |
     Discrim |   2.014221    .112486    17.91   0.000     1.793753     2.23469
        Diff |
         =1  |   .3911525   .0318698    12.27   0.000     .3286887    .4536162
-------------+----------------------------------------------------------------
ft6          |
     Discrim |   2.267844   .1059075    21.41   0.000      2.06027    2.475419
        Diff |
       >= 1  |  -1.086657     .04071   -26.69   0.000    -1.166447   -1.006867
       >= 2  |  -.9268428   .0381401   -24.30   0.000    -1.001596   -.8520896
       >= 3  |  -.7016476   .0349022   -20.10   0.000    -.7700547   -.6332406
       >= 4  |  -.4866712   .0325674   -14.94   0.000     -.550502   -.4228403
       >= 5  |  -.3082986   .0311084    -9.91   0.000    -.3692699   -.2473273
       >= 6  |   .1763741   .0295792     5.96   0.000        .1184    .2343483
       >= 7  |   .3869222   .0302634    12.79   0.000     .3276071    .4462374
       >= 8  |   .6974215   .0329238    21.18   0.000      .632892    .7619509
       >= 9  |   1.190076   .0413708    28.77   0.000     1.108991    1.271161
        =10  |   1.467087   .0487822    30.07   0.000     1.371476    1.562698
------------------------------------------------------------------------------

. predict thetapt, latent
(option ebmeans assumed)
(using 30 quadrature points)

. polychoric thetapt pref_pt iden_pt think_pt symp_pt, pw verbose 

Variables :  thetapt pref_pt
Type :       polyserial
Rho        = .93861736
S.e.       = .01148858

Variables :  thetapt iden_pt
Type :       polyserial
Rho        = .5317675
S.e.       = .02188712

Variables :  thetapt think_pt
Type :       polyserial
Rho        = .64531946
S.e.       = .0176608

Variables :  thetapt symp_pt
Type :       polyserial
Rho        = .67848259
S.e.       = .00656494
no observations

Variables :  pref_pt iden_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt symp_pt
Type :       polychoric
Rho        = .9187894
S.e.       = .01423103
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt symp_pt
Type :       polychoric
Rho        = .87293878
S.e.       = .02002008
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt symp_pt
Type :       polychoric
Rho        = .86335961
S.e.       = .02359722
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. drop thetapt

. irt grm pref_pt iden_pt symp_pt vote_pt ft6 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmethod(mcaghermite) 

Fitting fixed-effects model:

Iteration 0:   log likelihood = -11962.239  
Iteration 1:   log likelihood = -11962.239  

Fitting full model:

Iteration 0:   log pseudolikelihood = -11291.743  
Iteration 1:   log pseudolikelihood =  -10629.84  
Iteration 2:   log pseudolikelihood = -10600.776  
Iteration 3:   log pseudolikelihood = -10571.387  
Iteration 4:   log pseudolikelihood = -10569.272  
Iteration 5:   log pseudolikelihood = -10569.117  
Iteration 6:   log pseudolikelihood = -10569.124  
Iteration 7:   log pseudolikelihood = -10569.122  
Iteration 8:   log pseudolikelihood = -10569.122  

Graded response model                           Number of obs     =      4,324
Log pseudolikelihood = -10569.122
                         (Std. Err. adjusted for 3,118 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_pt      |
     Discrim |   8.479673   1.664384     5.09   0.000      5.21754    11.74181
        Diff |
         =1  |   .7249718   .0339671    21.34   0.000     .6583976     .791546
-------------+----------------------------------------------------------------
iden_pt      |
     Discrim |   5.326319   1.126234     4.73   0.000     3.118941    7.533698
        Diff |
         =1  |   1.002315   .0420897    23.81   0.000     .9198205    1.084809
-------------+----------------------------------------------------------------
symp_pt      |
     Discrim |   3.828023   .2991273    12.80   0.000     3.241745    4.414302
        Diff |
         =1  |    1.24073    .035638    34.81   0.000     1.170881    1.310579
-------------+----------------------------------------------------------------
vote_pt      |
     Discrim |    1.98758   .1095862    18.14   0.000     1.772795    2.202365
        Diff |
         =1  |   .3918868   .0320437    12.23   0.000     .3290823    .4546914
-------------+----------------------------------------------------------------
ft6          |
     Discrim |   2.257349   .1034435    21.82   0.000     2.054603    2.460094
        Diff |
       >= 1  |  -1.083224   .0405001   -26.75   0.000    -1.162603   -1.003846
       >= 2  |  -.9234268   .0379647   -24.32   0.000    -.9978362   -.8490174
       >= 3  |  -.6980731   .0347851   -20.07   0.000    -.7662506   -.6298956
       >= 4  |  -.4826996   .0324812   -14.86   0.000    -.5463616   -.4190376
       >= 5  |  -.3037631    .031021    -9.79   0.000    -.3645632   -.2429631
       >= 6  |   .1821146    .029579     6.16   0.000     .1241409    .2400884
       >= 7  |    .393021   .0302754    12.98   0.000     .3336823    .4523597
       >= 8  |   .7045659   .0329547    21.38   0.000     .6399758     .769156
       >= 9  |   1.201028   .0414032    29.01   0.000      1.11988    1.282177
        =10  |   1.481887   .0490509    30.21   0.000     1.385749    1.578026
------------------------------------------------------------------------------

. predict thetapt, latent
(option ebmeans assumed)
(using 30 quadrature points)

. polychoric thetapt pref_pt iden_pt think_pt symp_pt, pw verbose 

Variables :  thetapt pref_pt
Type :       polyserial
Rho        = .93821151
S.e.       = .01147878

Variables :  thetapt iden_pt
Type :       polyserial
Rho        = .76197387
S.e.       = .01758887

Variables :  thetapt think_pt
Type :       polyserial
Rho        = .48868432
S.e.       = .02296001

Variables :  thetapt symp_pt
Type :       polyserial
Rho        = .67169927
S.e.       = .00645123
no observations

Variables :  pref_pt iden_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt symp_pt
Type :       polychoric
Rho        = .9187894
S.e.       = .01423103
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt symp_pt
Type :       polychoric
Rho        = .87293879
S.e.       = .02002008
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt symp_pt
Type :       polychoric
Rho        = .86335961
S.e.       = .02359722
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. drop thetapt

. irt grm pref_pt iden_pt think_pt vote_pt ft6 if wave==1 | wave==6, intpoints(30) difficult vce(cluster idtelefone) intmethod(mcaghermite) 

Fitting fixed-effects model:

Iteration 0:   log likelihood = -10768.671  
Iteration 1:   log likelihood = -10768.671  

Fitting full model:

Iteration 0:   log pseudolikelihood = -10341.086  
Iteration 1:   log pseudolikelihood =  -9884.646  
Iteration 2:   log pseudolikelihood = -9877.8088  
Iteration 3:   log pseudolikelihood =  -9862.672  
Iteration 4:   log pseudolikelihood =  -9862.025  
Iteration 5:   log pseudolikelihood = -9862.0207  
Iteration 6:   log pseudolikelihood = -9862.0207  

Graded response model                           Number of obs     =      4,301
Log pseudolikelihood = -9862.0207
                         (Std. Err. adjusted for 3,106 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_pt      |
     Discrim |   3.739714    .537987     6.95   0.000     2.685279    4.794149
        Diff |
         =1  |   .7485864   .0441255    16.96   0.000      .662102    .8350709
-------------+----------------------------------------------------------------
iden_pt      |
     Discrim |   2.880531   .4014107     7.18   0.000     2.093781    3.667282
        Diff |
         =1  |   1.110986   .0619802    17.92   0.000     .9895073    1.232465
-------------+----------------------------------------------------------------
think_pt     |
     Discrim |   2.502805    .329682     7.59   0.000      1.85664     3.14897
        Diff |
         =1  |   1.454717   .0884896    16.44   0.000     1.281281    1.628154
-------------+----------------------------------------------------------------
vote_pt      |
     Discrim |   2.153929    .136134    15.82   0.000     1.887111    2.420747
        Diff |
         =1  |   .3739413   .0306224    12.21   0.000     .3139225    .4339601
-------------+----------------------------------------------------------------
ft6          |
     Discrim |   2.812167   .1931211    14.56   0.000     2.433656    3.190677
        Diff |
       >= 1  |  -1.014811   .0382088   -26.56   0.000    -1.089699   -.9399228
       >= 2  |  -.8677032   .0356167   -24.36   0.000    -.9375106   -.7978957
       >= 3  |  -.6597415   .0323971   -20.36   0.000    -.7232387   -.5962444
       >= 4  |  -.4603715   .0301056   -15.29   0.000    -.5193775   -.4013655
       >= 5  |  -.2944169   .0286817   -10.26   0.000     -.350632   -.2382018
       >= 6  |    .158887   .0273348     5.81   0.000     .1053119    .2124622
       >= 7  |   .3567991   .0282426    12.63   0.000     .3014447    .4121536
       >= 8  |   .6496633    .031284    20.77   0.000     .5883478    .7109789
       >= 9  |   1.118402   .0395651    28.27   0.000     1.040856    1.195948
        =10  |   1.385838   .0467048    29.67   0.000     1.294298    1.477378
------------------------------------------------------------------------------

. predict thetapt, latent
(option ebmeans assumed)
(using 30 quadrature points)

. polychoric thetapt pref_pt iden_pt think_pt symp_pt, pw verbose 

Variables :  thetapt pref_pt
Type :       polyserial
Rho        = .78709785
S.e.       = .01532202

Variables :  thetapt iden_pt
Type :       polyserial
Rho        = .64042448
S.e.       = .01533743

Variables :  thetapt think_pt
Type :       polyserial
Rho        = .54932187
S.e.       = .01830735

Variables :  thetapt symp_pt
Type :       polyserial
Rho        = .39711127
S.e.       = .00934712
no observations

Variables :  pref_pt iden_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_pt symp_pt
Type :       polychoric
Rho        = .9187894
S.e.       = .01423103
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_pt think_pt
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_pt symp_pt
Type :       polychoric
Rho        = .87293878
S.e.       = .02002008
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_pt symp_pt
Type :       polychoric
Rho        = .86335961
S.e.       = .02359722
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. 
.         *2: PSDB Results
. capture drop thetapsdb

. irt grm iden_psdb think_psdb symp_psdb vote_psdb ft7 if wave==1 | wave==6, intpoints(60) difficult vce(cluster idtelefone) intmethod(mcaghermite) from(star
> t, skip)

Iteration 0:   log pseudolikelihood = -7569.2581  
Iteration 1:   log pseudolikelihood = -7569.0365  
Iteration 2:   log pseudolikelihood = -7569.0321  
Iteration 3:   log pseudolikelihood = -7569.0321  

Graded response model                           Number of obs     =      4,322
Log pseudolikelihood = -7569.0321
                         (Std. Err. adjusted for 3,116 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
iden_psdb    |
     Discrim |    3.76565   .7495835     5.02   0.000     2.296494    5.234807
        Diff |
         =1  |   1.933923   .0969171    19.95   0.000     1.743969    2.123877
-------------+----------------------------------------------------------------
think_psdb   |
     Discrim |   3.426207    .694121     4.94   0.000     2.065755    4.786659
        Diff |
         =1  |   2.156022   .1320206    16.33   0.000     1.897266    2.414778
-------------+----------------------------------------------------------------
symp_psdb    |
     Discrim |   8.880557   3.914769     2.27   0.023     1.207751    16.55336
        Diff |
         =1  |   1.936301   .0821726    23.56   0.000     1.775246    2.097356
-------------+----------------------------------------------------------------
vote_psdb    |
     Discrim |    1.08436   .0983307    11.03   0.000     .8916354    1.277085
        Diff |
         =1  |   1.899045   .1329954    14.28   0.000     1.638379    2.159711
-------------+----------------------------------------------------------------
ft7          |
     Discrim |   1.329304   .1078165    12.33   0.000     1.117987     1.54062
        Diff |
       >= 1  |   -1.43838   .0875595   -16.43   0.000    -1.609994   -1.266767
       >= 2  |  -1.235159   .0780766   -15.82   0.000    -1.388187   -1.082132
       >= 3  |  -.8492625   .0617631   -13.75   0.000    -.9703158   -.7282091
       >= 4  |  -.4673943   .0486329    -9.61   0.000     -.562713   -.3720755
       >= 5  |  -.0711117   .0413249    -1.72   0.085     -.152107    .0098835
       >= 6  |   .7978913   .0563004    14.17   0.000     .6875446    .9082381
       >= 7  |   1.308594    .077527    16.88   0.000     1.156644    1.460544
       >= 8  |   1.879908   .1076406    17.46   0.000     1.668936     2.09088
       >= 9  |   2.753368   .1665775    16.53   0.000     2.426882    3.079854
        =10  |   3.147469   .1961516    16.05   0.000     2.763019    3.531919
------------------------------------------------------------------------------

. predict thetapsdb, latent
(option ebmeans assumed)
(using 60 quadrature points)

. polychoric thetapsdb pref_psdb iden_psdb think_psdb symp_psdb, pw verbose 

Variables :  thetapsdb pref_psdb
Type :       polyserial
Rho        = .3919923
S.e.       = .02084158

Variables :  thetapsdb iden_psdb
Type :       polyserial
Rho        = .42185212
S.e.       = .02495549

Variables :  thetapsdb think_psdb
Type :       polyserial
Rho        = .3537251
S.e.       = .02494017

Variables :  thetapsdb symp_psdb
Type :       polyserial
Rho        = .44842561
S.e.       = .00704181
no observations

Variables :  pref_psdb iden_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb symp_psdb
Type :       polychoric
Rho        = .96150611
S.e.       = .01442379
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb symp_psdb
Type :       polychoric
Rho        = .89534667
S.e.       = .03203524
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_psdb symp_psdb
Type :       polychoric
Rho        = .84448304
S.e.       = .04660182
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. drop thetapsdb

. irt grm pref_psdb think_psdb symp_psdb vote_psdb ft7 if wave==1 | wave==6, intpoints(60) difficult vce(cluster idtelefone) intmethod(mcaghermite) from(star
> t, skip)

Iteration 0:   log pseudolikelihood = -7572.8438  
Iteration 1:   log pseudolikelihood = -7572.7298  
Iteration 2:   log pseudolikelihood =  -7572.724  
Iteration 3:   log pseudolikelihood = -7572.7234  

Graded response model                           Number of obs     =      4,319
Log pseudolikelihood = -7572.7234
                         (Std. Err. adjusted for 3,115 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_psdb    |
     Discrim |   7.069785   2.442164     2.89   0.004     2.283231    11.85634
        Diff |
         =1  |   1.692514   .0452564    37.40   0.000     1.603813    1.781214
-------------+----------------------------------------------------------------
think_psdb   |
     Discrim |   3.392482   .6907894     4.91   0.000      2.03856    4.746405
        Diff |
         =1  |   2.160106   .1321648    16.34   0.000     1.901068    2.419144
-------------+----------------------------------------------------------------
symp_psdb    |
     Discrim |   9.519356   6.485799     1.47   0.142    -3.192576    22.23129
        Diff |
         =1  |   1.929755   .1047801    18.42   0.000      1.72439     2.13512
-------------+----------------------------------------------------------------
vote_psdb    |
     Discrim |   1.056725   .0972873    10.86   0.000     .8660455    1.247405
        Diff |
         =1  |    1.93342   .1403589    13.77   0.000     1.658322    2.208519
-------------+----------------------------------------------------------------
ft7          |
     Discrim |   1.343373   .1039582    12.92   0.000     1.139619    1.547128
        Diff |
       >= 1  |  -1.429631   .0838721   -17.05   0.000    -1.594018   -1.265245
       >= 2  |  -1.227524   .0749613   -16.38   0.000    -1.374445   -1.080603
       >= 3  |   -.843883   .0595639   -14.17   0.000    -.9606261     -.72714
       >= 4  |  -.4636627   .0477951    -9.70   0.000    -.5573394   -.3699861
       >= 5  |  -.0682317   .0412936    -1.65   0.098    -.1491656    .0127022
       >= 6  |   .7977118   .0556076    14.35   0.000     .6887229    .9067007
       >= 7  |   1.304911     .07574    17.23   0.000     1.156463    1.453358
       >= 8  |    1.87082   .1045068    17.90   0.000      1.66599    2.075649
       >= 9  |    2.73374    .160185    17.07   0.000     2.419783    3.047697
        =10  |   3.123761   .1915945    16.30   0.000     2.748243     3.49928
------------------------------------------------------------------------------

. predict thetapsdb, latent
(option ebmeans assumed)
(using 60 quadrature points)

. polychoric thetapsdb pref_psdb iden_psdb think_psdb symp_psdb, pw verbose 

Variables :  thetapsdb pref_psdb
Type :       polyserial
Rho        = .49387108
S.e.       = .02171424

Variables :  thetapsdb iden_psdb
Type :       polyserial
Rho        = .32890985
S.e.       = .02559452

Variables :  thetapsdb think_psdb
Type :       polyserial
Rho        = .35433539
S.e.       = .02518926

Variables :  thetapsdb symp_psdb
Type :       polyserial
Rho        = .44235513
S.e.       = .00756743
no observations

Variables :  pref_psdb iden_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb symp_psdb
Type :       polychoric
Rho        = .96150612
S.e.       = .01442379
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb symp_psdb
Type :       polychoric
Rho        = .89534667
S.e.       = .03203524
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_psdb symp_psdb
Type :       polychoric
Rho        = .84448304
S.e.       = .04660182
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. drop thetapsdb

. irt grm pref_psdb iden_psdb think_psdb vote_psdb ft7 if wave==1 | wave==6, intpoints(60) difficult vce(cluster idtelefone) intmethod(mcaghermite) from(star
> t, skip)

Iteration 0:   log pseudolikelihood =  -7444.761  (not concave)
Iteration 1:   log pseudolikelihood = -7443.9169  
Iteration 2:   log pseudolikelihood =  -7443.662  
Iteration 3:   log pseudolikelihood = -7443.5686  
Iteration 4:   log pseudolikelihood = -7443.5614  
Iteration 5:   log pseudolikelihood = -7443.5613  

Graded response model                           Number of obs     =      4,287
Log pseudolikelihood = -7443.5613
                         (Std. Err. adjusted for 3,095 clusters in idtelefone)
------------------------------------------------------------------------------
             |               Robust
             |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
pref_psdb    |
     Discrim |   4.669701   2.407091     1.94   0.052    -.0481108    9.387512
        Diff |
         =1  |   1.731581    .140348    12.34   0.000     1.456504    2.006658
-------------+----------------------------------------------------------------
iden_psdb    |
     Discrim |   3.092123   1.120252     2.76   0.006     .8964689    5.287777
        Diff |
         =1  |   2.031518   .2109439     9.63   0.000     1.618076    2.444961
-------------+----------------------------------------------------------------
think_psdb   |
     Discrim |    5.27654   3.883548     1.36   0.174    -2.335074    12.88815
        Diff |
         =1  |   2.034801    .190586    10.68   0.000      1.66126    2.408343
-------------+----------------------------------------------------------------
vote_psdb    |
     Discrim |    1.10758   .1133202     9.77   0.000     .8854762    1.329683
        Diff |
         =1  |   1.869656   .1419573    13.17   0.000     1.591425    2.147887
-------------+----------------------------------------------------------------
ft7          |
     Discrim |   1.300526   .1345551     9.67   0.000     1.036803    1.564249
        Diff |
       >= 1  |  -1.454427   .1062522   -13.69   0.000    -1.662677   -1.246176
       >= 2  |  -1.248112    .093576   -13.34   0.000    -1.431517   -1.064706
       >= 3  |  -.8564962   .0708038   -12.10   0.000    -.9952691   -.7177233
       >= 4  |  -.4690351   .0521616    -8.99   0.000    -.5712699   -.3668003
       >= 5  |  -.0668067   .0417062    -1.60   0.109    -.1485494     .014936
       >= 6  |    .814409   .0656839    12.40   0.000      .685671     .943147
       >= 7  |   1.329741   .0949595    14.00   0.000     1.143624    1.515858
       >= 8  |    1.90612   .1334441    14.28   0.000     1.644574    2.167665
       >= 9  |   2.794232    .205011    13.63   0.000     2.392418    3.196046
        =10  |   3.198118   .2436132    13.13   0.000     2.720644    3.675591
------------------------------------------------------------------------------

. predict thetapsdb, latent
(option ebmeans assumed)
(using 60 quadrature points)

. polychoric thetapsdb pref_psdb iden_psdb think_psdb symp_psdb, pw verbose 

Variables :  thetapsdb pref_psdb
Type :       polyserial
Rho        = .45654467
S.e.       = .01239888

Variables :  thetapsdb iden_psdb
Type :       polyserial
Rho        = .42523633
S.e.       = .0155693

Variables :  thetapsdb think_psdb
Type :       polyserial
Rho        = .3689515
S.e.       = .01529943

Variables :  thetapsdb symp_psdb
Type :       polyserial
Rho        = .26674544
S.e.       = .01113687
no observations

Variables :  pref_psdb iden_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
insufficient observations

Variables :  pref_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  pref_psdb symp_psdb
Type :       polychoric
Rho        = .96150623
S.e.       = .01442374
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .
could not calculate numerical derivatives
missing values encountered
could not calculate numerical derivatives
missing values encountered
no variables defined

Variables :  iden_psdb think_psdb
Type :       polychoric
Rho        = .
S.e.       = .
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  iden_psdb symp_psdb
Type :       polychoric
Rho        = .89534667
S.e.       = .03203524
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

Variables :  think_psdb symp_psdb
Type :       polychoric
Rho        = .84448304
S.e.       = .04660182
Goodness of fit tests:
Pearson G2 = ., Prob( >chi2(.)) = .
LR X2      = ., Prob( >chi2(.)) = .

. 
. log close
      name:  <unnamed>
       log:  C:\Users\Andy baker\OneDrive - University of Colorado at Boulder Office 365\My Documents\Research\Paper PID Experiment\Brazil\JOP IRT Brazil Res
> ults.log
  log type:  text
 closed on:   7 Mar 2018, 01:35:18
-------------------------------------------------------------------------------------------------------------------------------------------------------------
